Conductivity of graphene: How to distinguish between samples with short and long 

range scatterers 
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Applying a quasiclassical equation to carriers in graphene we found a way how to distinguish 
between samples with the domination of short and long range scatterers from the conductivity 
measurements. The model proposed explains recent transport experiments with chemically doped 
as well as suspended graphene. 



INTRODUCTION 

To understand the spectacular transport properties of 
single graphite layers [l|, [2| one has to know which kind 
of scatterers dominates in a given sample. Indeed, the- 
ory predicts the electrical conductivity to be strongly de- 
pendent on the particular type of the scatterers being 
present in the system Q • We demonstrate that the diver- 
sity of the conductivity behaviour observed in graphene 
0, 0, @, 0, 0] can be naturally described employing the 
concentrations of long and short range scatterers as pa- 
rameters. According to our microscopic model, the con- 
ductivity measurements upon potassium doping [^| ob- 
viously suggest a strong domination of the long range 
scattering potential at all possible carrier concentrations, 
whereas the nitrogen dioxide doping [f| induces charged 
impurities dominating the scattering mechanisms near 
the conductivity minimum only. In addition, we shall 
see that zero chemical doping gives rise to the very sharp 
minimum conductivity dip observed in Refs. Our 
semiclassical approach fully incorporates the chiral na- 
ture of electronic states in graphene which is crucial for 
a correct description of transport properties. 

There are three most important regimes to consider: 
(I) — Long range scattering potential limit. The long 
range potential of charged impurities strongly dominates 
the scattering processes and governs both the carrier mo- 
bility and the conductivity minimum, as shown in fig. [1] 
and observed in Ref. Here, the width of the minimum 
conductivity region is rather broad, with a well defined 
plateau, and the conductivity minimum itself is around 
the universal value 4e 2 //i. Upon further chemical dop- 
ing the carrier mobility decreases, the conductivity as 
a function of carrier concentration above its minimum 
becomes more linear, whereas the residual conductivity 
defined by the linear fit keeps constant near 2e 2 /h. (II) 
- Both long and short range scattering potentials are 
in play. At somewhat lower charged impurity concen- 
tration, the long range scatterers do not play an essen- 
tial role in the conductivity behaviour far away from the 
conductivity minimum. Here, the short range scatter- 
ers (i. e. nano-sized ripples [9( and other imperfections 
[loj ]) govern the carrier mobility which is defined by the 
linear fit as shown fig. [2] and, thus, turns out to be in- 
sensitive to the chemical doping, as observed in Ref. 



The minimum conductivity dip is rather sharp, i. e. no 
plateau is unambiguously visible. The conductivity min- 
imum is still around 4e 2 //i, whereas the residual con- 
ductivity can acquire a wide range of values depending 
on the scattering parameters. (Ill) — Short range scat- 
tering potential limit. If the concentration of the long 
range scatterers vanishes then the minimum conductiv- 
ity dip becomes verysharp as it is shown in fig. [3] and 
observed in Refs. Moreover, at zero concentration 

of long range scatterers the conductivity minimum be- 
comes indistinguishable from its residual value which is 
not universal and can reach the values up to tens of e 2 /h 
as observed in Ref. Thus, in this latter case both 
the minimum conductivity dip and carrier mobility are 
governed by the short range scatterers. 

To describe these three cases theoretically we choose 
a quasiclassical Boltzmann approach since it allows us 
to get transparent analytical formulas for the electrical 
conductivity in a broad range of parameters. The most 



of the previous attempts [llj, [12|, [13|, [1J, [15| based on 



the Boltzmann equation discard the chirality of low en- 
ergy excitations, and we argue that a proper model must 
necessary contain both peculiar properties of carriers in 
graphene, namely, the linear spectrum of low energy exci- 
tations and its chiral nature. The quasiclassical approach 
looks at the first sight inapplicable to investigate the con- 
ductivity near its minimum since it is expected to fail 
at low Fermi energies Ep as soon as Eft(Ef) becomes 
comparable with h. Note, however, that the relaxation 
time t for short range scatterers diverges 3 at Ep = 
that leads to a constant product Eft(Ef) larger than H. 
Moreover, as we shall see below, for the long range impu- 
rity potentials the conductivity approaches its minimum 
away from the Dirac point which makes the quasiclassical 
model applicable again. 

We start from the Dirac Hamiltonian Hq — hvo(a x k x + 
o~yky) describing the low energy excitations in the ir- 
system of graphene around the K corner of the first Bril- 
louin zone. [Here vq w 10 6 ms _1 is the effective "speed of 
light" , o~ x ,y are the Pauli matrices describing the pseudo- 
(or sublattice-)spin, and k is the two-component particle 
momentum.] To simplify the model we adopt the valley 
and spin degeneracy of the electron states, see Appendix 
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for details. The eigenstates have the form 
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where tan# = k y /k x , and the energy spectrum reads 
Ek± — ±Hvok. The spinors |T]) contain momentum- 
dependent phase factor e lS which entangles the momen- 
tum and the sublattice degree of freedom. This chirality 
of electronic states plays a crucial role in our description. 

To consider both short and long range scatterers on 
the one hand and to avoid having too many fitting pa- 
rameters on the other, we choose two opposite types of 
scattering potential, namely, Uq(t) = lira U(r) (with 

Q^R 

R being the radius of the short range scatterers), and 
t/i(r) = lim U(r), where U(r) is the generic one given 

g — >oc 

by U(r) — (u /r)c~ r / e . Assuming that the short and 
long range scatterers have the concentrations No and Ni 
respectively, it is easy to derive the corresponding relax- 
ation times Q given by 
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where no,i = ttNo.i(2uo / hvo) 2 are the renormalised con- 
centrations being only the two fitting parameters. The 
characteristic size R of ripples and other imperfections is 
assumed to be near a few nm. 

To incorporate the chiral nature of Dirac fermions into 
the stationary quasiclassical kinetic equation linear in the 
homogeneous electric field E we have to start from its 
general form 



1 / „df 



H J(k) = ![/], 
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where / is the density matrix, is the collision in- 
tegral. In eq. the sublattice degree of freedoom is 
assumed to be a quantum number, although the particle 
momentum k is still considered quasiclassically. To deal 
with such an equation it is natural to rewrite it in the 
basis of the eigenvectors (TTJ). The distribution function 
/ represents then a 2 x 2 matrix, and its off-diagonal el- 
ements do not drop out thanks to the chiral nature of 
the eigenstates (J]) . This important modification of the 
Boltzmann equation allows us to deal with the essentially 
quantum effect suggesting that a particle could not only 
be in one of the states v E'k+ or Wk- but in an arbitrary su- 
perposition of them. Alternatively one might think about 
the off-diagonal elements as a manifestation of the pseu- 
dospin precession 0] associated with the chirality index 

K. 

To solve the kinetic equation with respect to / we as- 
sume the linear response regime (small E). All techni- 
cal details concerning the kinetic equation, its solution, 



and derivation of conductivity expressions are given in 
the Appendix section. To write down the final zero- 
temperature conductivity formula it is convenient to in- 
troduce the following parameters: 70 = TrR 2 n and 
7„ = ■nR 2 n, where n — kp/n is the carrier concentra- 
tion with kp = Ep/hvo being the Fermi wave vector. 
Then, the conductivity reads 
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The sharp peak at n = in eq. () due to the term 
cx 2ni/n should be seen as an artifact since our quasi- 
classical approach certainly fails in this limit. According 
to the experimental data there is normally a plateau in- 
stead of that peak. We define the residual conductivity 
<r'°" g from the linear fit for the conductivity measure- 
ments in the limit of long range scatterers domination, 
as explained in Ref. As follows from the measure- 
ments [4], crj;™ s turns out to be a constant nearly equal 
to 2e 2 /h. The residual conductivity cr 1 ^ 6 is the only pa- 
rameter that could not be derived within our quasiclassi- 
cal approach because of the obvious restriction discussed 
before. But it contributes to cr(n) in a very simple way: 
<t'°" s = 2e 2 /h just shifts all the cr(n) curves into the 
region of higher conductivities and does not have an in- 
fluence on the quasiclassical physics above cr'°" g . Note, 
however, that the residual conductivity is not a constant 
as soon as the long range scatterers do not dominate any 
more. Such residual conductivity Crcs 01 * > °r°s S can be 
derived within our approach since it remains valid for 
short range scatterers even at zero carrier concentration. 
Note, that eq. () is true for arbitrary values of no,i and R, 
thus, we can utilise this formula to distinguish the trans- 
port regimes with the domination of short or long range 
scatterers. To proceed we set R = lnm, thus 7„ <C 1 
for all reasonable carrier concentrations, and, as conse- 
quence, 70 <C 1 as well. 



(I) — LONG RANGE SCATTERING POTENTIAL 
LIMIT 



Long range scatterers dominate if ri\ > no- Here, eq. () 
can be expanded in terms of 7„ and 7 n ^o/ n i resulting in 
the following expression 
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FIG. 1: Conductivity vs. carrier concentration: Domina- 
tion of the long range scatterers with the renormalised con- 
centration ni [defined below eqs. (J2j{3j)] which increases for 
each curve from 0.1 (upper curve) to 0.5 [xl0 12 cm -2 ] with 
0.1 ■ 10 12 cm -2 step. The renormalised concentration of the 
short range scatterers no [defined below eqs. ([2]-[3])] keeps con- 
stant at 0.1 TO 12 cm -2 . Dotted lines are the linear approxima- 
tions given by cr long = e^ long n + cr^ lg where ^ long = e/{2hm) 
(with h — 2nh being the Planck constant) is the carrier mo- 
bility obviously governed by the long range scatterers. The 
conductivity minimum a m i n is close to 4e 2 //i whereas the 
residual conductivity a TCB is equal to 2e 2 /h and turns out 
to be independent of m. The sharp peak at n = due to 
the term oc 2ni/n in eq. () has been removed as an artifact. 
This transport regime corresponds to the graphene samples 
strongly doped by potassium Q. 




Carrier concentration [10 cm ] 



FIG. 2: Conductivity vs. carrier concentration: Both long 
and short range scattering potentials contribute essentially. 
The renormalised concentration of short range scatterers no 
is assumed to be a constant 10 12 cm -2 , whereas the concentra- 
tion of long range scatterers m increases for each curve from 
zero (upper curve) to 0.25 ■ 10 12 cm~ 2 ] with 0.05 ■ 10 12 cm" 2 
step. Dotted lines are the linear approximations given by 
cr short = eu short n + cr* ort with the carrier mobility u short = 
e/(hno) governed by the short range scatterers and <7res° rt de- 
fined by fit. Though the conductivity minimum is still near 
4e 2 /h and does not depend strongly on the impurity concen- 
tration ni , the residual conductivity depends on both concen- 
trations for short and long range scatterers, see the main text. 
This transport regime corresponds to the graphene samples 
doped by nitrogen dioxide [5J. 
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The function er long (n) has a minimum near n m in = 2ni, 
rii 0. At this concentration the term proportional to 70 
is much smaller then the leading one, and the conductiv- 
ity reaches its minimum close to o\°" s = er!°" s + 2e 2 /h as 

" IYIIII ICO ' 

it is seen in fig. [1] Since erj™ g = 2e 2 /h the conductivity 
minimum indeed acquires the typical value 4e 2 / h in quite 
perfect samples (i. e. without ripples, neutral impurities 
etc.) when 70 — > 0. The corrections of the order of 0(70) 
and 0(70 In 7„) give all together a small negative contri- 
bution to the conductivity minimum. Most importantly, 
the conductivity minimum does not depend on m, i. e. 
it should be insensitive to doping, whereas the minimum 
carrier concentration rimin = 2rti increases linearly with 
chemical doping and in that way gives rise to the width 
of the minimum conductivity plateau, as observed in 0] • 
Above n = n m i n (i. e. n 3> 2ni) one can linearise eq. (O 
and obtain cr long = e/i lon Sn + crJ.™ g with /j} oas = e/{2hn l ) 
being the carrier mobility. The deviations from linear de- 
pendency cr long (n) are described by the term oc —n 2 /2n\ 
in eq. {fij and can be seen in fig. [T] as well. All these 
peculiar features indicate the domination of long range 
scatterers that was observed in Q. 



(II) — BOTH LONG AND SHORT RANGE 
SCATTERING POTENTIALS ARE IN PLAY 



At lower charged impurity concentrations the mini- 
mum conductivity dip becomes sharper, the plateau van- 
ishes, and the carrier mobility is governed by the short 
range scatterers, see fig. [2j For large enough carrier con- 
centrations so that no7n 3> n\ the conductivity can be 
written as 
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We emphasise that eq. © is just an estimate for the 
conductivity far away from its minimum, i. e. for the 
carrier concentrations n 3> ^1/70- Eq- © allows us to 
define the carrier mobility as ^t short — e/(hno) and resid- 

ual conductivity ^ = T 4 (l " ^ 
7at being j n taken at n = N large enough to satisfy 
N » ni/70. Though the estimation for CT^ s ort is very 
rough [since eq. ijHJ) is an estimate as well], it shows that 
the residual conductivity is not a constant anymore as 
long as the short range scatterers dominate. Thus, the 
residual conductivity depends on both short and long 
range scatterers, whereas the conductivity minimum is 
governed by the long range scatterers alone and remains 
near the typical value Ae 2 /h, as one can observe in fig. [2] 
The latter is due to the fact that n <C t at vanishing 
carrier density even if «o > nj. On the other hand, the 
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FIG. 3: Conductivity vs. carrier concentration: Domination 
of the short range scatterers with the renormalised concentra- 
tion no = 10 12 cm -2 . The renormalised concentration of long 
range scatterers m increases for each curve from 5 • 10 9 cm -2 
(upper curve) to 2.5 • 10 10 cm" 2 with 5 ■ 10 9 cm" 2 step. The 
dotted line corresponds to the case n\ = 0. The conductivity 
dip is very sharp and does not contain a plateau. According 
to the very recent report 0, the conductivity of suspended 
graphene after annealing indeed demonstrates very sharp dip 
at zero gate voltages. Comparison of the conductivity mea- 
surements for suspended and non-suspended graphene [g| is 
also consistent with our predictions. 

carrier mobility ^ short depends on no and, thus, is unaf- 
fected by the chemical doping. Such behaviour indicates 
the short range scatterers domination in the conductivity 
measurements upon nitrogen dioxide doping 

(Ill) — SHORT RANGE SCATTERING 
POTENTIAL LIMIT 

It is instructive to consider the special case n\ = in 
detail. Here the term oc 2n\/n does not dominate over 
the logarithmic one oc ln7„ at n — > 0, and, therefore, the 
latter must be retained. The conductivity then reads 

= T~ (w~ + — + 270 ln 7 „) + (9) 

Interrestingly, the conductivity minimum is uncertain, 
but since R is small (70, n <C 1) the conductivity at 
11 — ► is indeed indistinguishable from its residual value 
cr^° rt (ni = 0) = ^-^ + (i'°s s - This case is depicted in 
fig. [3] and corresponds to the chemically undoped @ and 
suspended @, H| graphene samples where the conductiv- 
ity demonstrates rather sharp dip with nonuniversal min- 
imum value strongly dependent on scattering parameters. 
It is noteworthy that similar logarithmic corrections to 
the conductivity minimum have been found in presence 
of electron-electron interactions [161 | l 71| and scattering 
on vacancies, cracks, or boundaries [l4j. 



Given the above theoretical results, let us now analyse 
recent experiments on graphene samples upon chemical 
doping [J, Q • We first discuss the case of potassium dop- 
ing [4j. Here, the width of the minimum conductivity 
region broadens (i. e. n m ; n increases) whereas both the 
minimum conductivity and mobility decrease, at least ini- 
tially. At certain stage, however, the conductivity min- 
imum becomes insensitive to further doping. These fea- 
tures are clearly understood within our model. Initially, 
the sample is very weakly doped, and the minimum con- 
ductivity is governed by the short range scatterers rather 
than charged impurities in accordance with eq. ([9]). At 
larger potassium concentrations the long range scatterers 
start to dominate, the minimum concentration n min oc ri\ 
increases whereas the conductivity minimum keeps con- 
stant despite further decreasing the mobility. This limit 
corresponds to eq. (J7J). In Ref. graphene is dosed 
by nitrogen dioxide. Here, the width of the conductivity 
minimum region broadens similar to the previous case 
but it remains much more narrow than upon potassium 
doping. Moreover, the carrier mobility and conductiv- 
ity minimum are both insensitive to doping. According 
to eq. ([8]) such data suggest strong domination of the 
short range scatterers, at least at higher carrier concen- 
trations n 3> rai/70. The reason of such a strong suppres- 
sion of induced charges might be the strong screening of 
Coulomb interactions in the samples examined. Such pre- 
cise description of charged impurities is certainly beyond 
the scope of this report, where U\ is not screened at all. 
We can predict, however, that the minimum conductivity 
phenomenon takes place even in free hanging graphene 
sheets, without any influence of the silicon dioxide sub- 
strate considered in [l5j |. 

Moreover, our model offers a solution to the urgent 
minimum conductivity problem. Previous measurements 
appear to be contradictory in so far as quite a represen- 
tative group of graphene samples [l[ exhibits the conduc- 
tivity minimum clustering around / Le 2 /h, whereas many 
other samples @ demonstrate very wide range of the 
minimum conductivity values, away from 4e 2 //i. Now 
it is clear that the universal conductivity minimum in 
the samples from [l[ is due to the long range scatterers 
(e. g. charged impurities) with a little influence of rip- 
ples and other short range imperfections. This results in 
the conductivity dip broaden enough to resolve its bot- 
tom at cr ~ Ae 2 /h. (The carrier mobility, however, is not 
necessarily governed by charged impurities.) In contrast, 
the conductivity minimum far away from 4e 2 jh betrays 
the lack of the charged impurities in a given sample that 
makes the conductivity dip sharper and its bottom dif- 
ficult to resolve, as it is seen from eq. ([9]). Thus, the 
minimum conductivity becomes indistinguishable from 
its residual value which, in general, is sensitive to the 
concentrations of both long and short range scatterers. 
It is therefore not surprising that the conductivity min- 
imum changes from sample to sample in this case. In 
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that way our model is able to explain the wide diversity 
of the minimum conductivity values measured in single 
layer graphenc. 

As an improtant further conclusion we would like to 
emphasise that it is the chiral nature of low energy excita- 
tions, rather than its linear spectrum, that is responsible 
for the minimum conductivity phenomenon in graphene. 
It is especially obvious for quite perfect graphenc sam- 
ples where long range scatterers (i. e. charged impuri- 
ties) dominate the scattering processes: The conductivity 
would be given just by the first term in eq. ([7]), and the 
conductivity minimum would certainly vanish as long as 
the chiral nature were discarded. Such a conclusion is 
consistent with the theoretical investigation of ballistic 
graphene [HI and strongly supported by the very ex- 
perimental fact that bilayer graphene, despite obviously 
non-linear dispersion law, also exhibits a minimum con- 
ductivity of the same order as single layer samples 
The conductivity measurements for chemically doped bi- 
layer samples are not accessible yet, they would certainly 
be an experimentum crucis on which the mechanism is 
responsible for the conductivity minimum in graphenc. 



APPENDIX — SOLUTION OF THE KINETIC 
EQUATION 

The distribution function /(k) for carriers in presence 
of the electric field E can be written as the sum of an 
equilibrium contribution f°(E kK ) (which is just the Fermi 
function) and a nonequilibrium part / x (k) which repre- 
sents a 2 x 2 matrix. To find / 1 (k) one has to solve 
eq. ([4]) written in the helicity basis To do that one 
has to bear in mind that the spinors are k-dependent (via 
the phase i8), thus, the kinetic equation differs from its 
standart analogue and reads 



degree of freedom couples electron and hole states via the 
off-diagonal terms Hq k in the Hamiltonian Ho that leads 
to the non-diagonal velocity operator and electron-hole 
coherent contributions in the conductivity. 

To find the collision integral in eq. (fTU|) we start from 
the commutator — i \U, f] , where the density matrix is 
substituted by the following expression describing its 
time evolution after each scattering event 



f -> /(* = 0) + lJ dte-i H ^[f(t = 0), U]e* Hot . (12) 



The integrals over t can be taken introducing a small 
parameter in the time dependent exponent 
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where P denotes the principle value. Then the collision 
integral reads 



= J 
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with Vy being the velocity matrix elements given by 
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i sin 8 — cos 8 
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(11) 

One might ask why spin and valley degrees of freedom 
are discarded in eq. (fTQ|) . The reason is that the total 
Hamiltonian H (which includes, in addition to the sub- 
lattice index k, the spin s and valley v indices) does not 
contain the crossing (off-diagonal) terms H ss and H vv , 
and, therefore spin and valley coherence does not mani- 
fest itself in the conductivity. In contrast, the sublattice 



i=0,l 



This form of the collision integral describes elastic scat- 
tering of chiral particles within Born approximation and 
is routinely applied to kinetics of chiral electrons with 

HE" 



spin-orbit interactions, see e.g. [191. 1201]. It is important 
that the intrinsicly inelastic terms [i. e. the principle 
value in eq. (fT3j) ] are discarded in the final expression for 
/[/]. These terms describe creation of electron- hole pairs 
and have been dubbed in 21J as zitterbewegung contri- 
butions. The latter does not change results qualitatively 
but essentially complicates the analytics as one can see 



comparing Refs. [2l| and [22j. To reach a compromise 



between transparency and completeness of our model we 
deal only with the leading (i. e. <5-functional) terms in 
the collision integral. 

The solution of eq. (|10[) can be written down as /(k) = 
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and /21 can be obtained from eqs. P^HTB|) just ex- 
changing the indices belong to E k and v accordingly. 



is the total momentum relaxation 



2 /ft 2 is the electron- hole inco- 
The off-diagonal elements f\ 2 and 



Here, = +t- 
time, and a(k) 
herence parameter 
/21 do not contribute to the current if and only if they 
become real at a ^> 1. 

The electrical current j = ej rfpTr [v(fc)/ 1 (k)] is 
easy to derive at low temperatures (T <C -Ef), and the 
conductivity can be written as a sum a = j/E + cr].™ 8 
with a constant crj.™ 8 being the residual conductivity in 
the limit of the long range scatterers domination defined 
by the linear fit from the experimental data [3] • The for- 
mula () follows directly from eqs. ((T5TTT6]l and is exact in 
this sense. It is only essential to integrate over k before 
the subsequent approximations relying on a small i?; we 
would obtain the logarithmic divergence in the conduc- 
tivity integral otherwise 21 , 23] . 
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